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Experimental and theoretical investigations of undulation patterns in high-pressure, inclined 
layer gas convection at a Prandtl number near unity are reported. Particular focus is given to the 
competition between the spatiotemporal chaotic state of undulation chaos and stationary patterns 
of ordered undulations. In experiments a competition and bistability between the two states is ob- 
served, with ordered undulations most prevalent at higher Rayleigh number. The spectral pattern 
entropy, spatial correlation lengths, and defect statistics are used to characterize the competing 
states. The experiments are complemented by a theoretical analysis of the Oberbeck-Boussinesq 
equations. The stability region of the ordered undulation as a function of their wavevectors and 
the Rayleigh number is obtained with Galerkin techniques. In addition, direct numerical sim- 
ulations are used to investigate the spatiotemporal dynamics. In the simulations both ordered 
undulations and undulation chaos were observed dependent on initial conditions. Experiment 
and theory are found to agree well. 



1. Introduction 

Many spatially extended pattern-forming systems show non-transient steady states in which 
spatial and temporal correlations fall off quickly; such states are often described as exhibiting 
"spatiotemporal chaos" (STC). Due to the coupling of spatial and temporal degrees of free- 
dom, STC is richer than the purely temporal chaos observed in many low-dimensional nonlinear 
systems. While spatiotemporally chaotic phenomena have received a great deal of attention in 
recent years, our understanding of STC is far from satisfactory (Cross and Hohenberg 1994; 
Gollub 1994; Egolf et al. 1998, 2000; Gollub and Cross 2000; Bodenschatz et al. 2000; Daniels 
and Bodenschatz 2002a). One important issue is to reveal universal mechanisms underlying STC 
phenomena in quite diverse systems. In this paper, we focus on striped patterns observed, for in- 
stance, in the shifting patterns of sand dunes, in cloud street formations, and in many biological 
systems (Cross and Hohenberg 1993). A well-investigated manifestation of STC in this setting 
is commonly described as defect turbulence, in which the creation, annihilation, and motion of 
topological point defects continually change the local wavenumber and orientation of the stripes 
(Ramazza et al. 1992; Rehberg et al. 1989; Bodenschatz et al. 2000; Daniels and Bodenschatz 
2002a; Young and Riecke 2003). Another feature is the bistability and competition between STC 
and well-ordered structures. Both may coexist in the same experimental snapshot, while in other 
cases the system switches between the two chaotically. While such a scenario seems to be com- 
mon in many pattern forming systems, only a few cases have been investigated in detail (Cakmur 
et al. 1997; Echebarria and Riecke 2000). 

Here, we study anisotropic stripe (roll) patterns in thermal convection of an inclined fluid layer 
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FIGURE 1 . Sample shadowgraph images of a convection cell inclined by an angle 7 = 30° and at Prandtl 
number P — 1.1 ± 0.04 with reduced driving e = 0.17 above onset. Uphill is left; the images have been 
Fourier-filtered to remove higher harmonics, (a) Undulation chaos (UC) and (b) ordered undulations (OU) 
(Daniels and Bodenschatz 2002a,b) 

(Daniels et al. 2000). In particular, we present a detailed experimental study of the bistability 
between ordered undulations (OU) of stripes and the defect turbulent state of undulation chaos 
(UC)) (see Fig. 1). The analysis is supported by numerical simulations (Brausch 2001) of the 
full Oberbeck-Boussinesq equations (OBE) where the OU and UC attractors were explored in 
a controlled manner. We characterize these two states in terms of the spectral pattern entropy, 
spatial correlation length, and defect density. These measures are found to be correlated with 
each other and are well suited to distinguish UC from OU. Along with these standard measures 
of STC, we determine the local wavevector of the patterns and describe the undulations of the 
stripes in terms of three characteristic amplitudes. These additional quantities allow for a refined 
discrimination of OU and UC and fit well into the theoretically determined stability island of 
undulations. 

2. Inclined Layer Convection 

Rayleigh-Benard convection (RBC, the thermal instability of a thin horizontal fluid layer 
heated from below) and its variants have been particularly fruitful in the investigation of pattern 
formation in extended systems (Bodenschatz et al. 2000). Experimentally, the use of compressed 
gases has allowed for the construction of shallow convection cells with large lateral extent and 
fast time scales (de Bruyn et al. 1996; Bodenschatz et al. 2000). In addition, the underlying fluid 
dynamical equations (Oberbeck-Boussinesq) are well-established and numerical techniques have 
been developed to simulate systems with many convection rolls (Pesch 1996; Bodenschatz et al. 
2000). 

One natural variation of RBC is inclined layer convection (ILC), in which the thin fluid layer is 
additionally tilted by an angle 7 with respect to horizontal, making the system anisotropic. This 
situation is common in nature as convective systems are often inclined with respect to gravity. As 
shown schematically in Fig. 2, the component of gravity parallel to the layer generates a cubic 
shear flow profile, Uo(z), upwards along the hotter (lower) plate and downward along the cooler 
(upper) plate. Therefore, ILC may also serve as a prototype system for convective systems in the 
presence of shear flow. 

As in RBC, the fluid becomes unstable to convection rolls above a critical temperature differ- 
ence, AT C . We use the reduced control parameter e = — 1 to measure the distance from the 
primary instability for fixed angle of inclination 7. For a Prandtl number P — vjn m 1 (with 



Ordered undulations and undulation chaos 



3 



A 




FIGURE 2. Schematic of a convection cell of thickness d inclined by an angle 7 subject to a temperature 
difference AT = T\ — T2 with associated coordinate system. Uo indicates the induced shearfiow (see 

text). 
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FIGURE 3. Phase diagram for inclined layer convection at P = 1.07, after (Daniels et al. 2000). Data 
points (and solid lines as guides to the eye) are observed boundaries between the different morphologies. 
Dashed line is Galerkin prediction for the instability of LR to OU. Gray arrow represents the region of data 
collection in this paper. 

kinematic viscosity v and thermal diffusivity n) both buoyancy and shear driven instabilities are 
observed, which evolve into numerous spatiotemporally chaotic states (Daniels et al. 2000). In 
Fig. 3, we have reproduced the part of the phase diagram relevant for the patterns examined in 
this paper. 

At small angles of inclination, buoyancy provides the primary instability and the convection 
rolls are aligned with the shear flow direction (longitudinal rolls, LR) (Clever 1973). Above a 
critical angle -f c2 (~ 78° for the parameter regime presented here), the primary instability is 
shear flow driven and the rolls align perpendicularly to the shear flow direction (transverse rolls, 
TR) (Hart 1971a,b; Clever and Busse 1977; Ruth 1980; Ruth et al. 1980). 

Over a range of intermediate angles (15° < 7 < 70°), transverse modes trigger a secondary 
bifurcation of LR to three dimensional undulation patterns (Clever and Busse 1977) slightly 
above onset (see the dashed line in Fig. 3 at e « 0.015.) A tertiary instability to a state of 
crawling rolls (at e m 0.3) limits the existence region of undulations from above. For a more 
detailed discussion of the phase diagram, see Daniels et al. (2000). 

Experimentally, we observed chaotic undulation patterns in the regime between the lower 
dashed and the upper solid destabilization line ("undulation chaos"), where undulating convec- 
tion rolls perpetually break and reconnect via moving point defects. Here, we restricted our in- 
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vestigations to a fixed angle of 7 = 30° for a range of e well inside the stability regime shown in 
Fig. 3. 

ILC can serve as a paradigm for a class of anisotropic pattern-forming systems such as liquid 
crystal convection (Kramer and Pesch 1995), Taylor-Couette flow (Tagg 1994), annular con- 
vection (Kurt et al. 2004), sand ripples (Blondeaux 1990; Hansen et al. 2001), or optical pattern 
formation (Ramazza et al. 1992) which also exhibit defect turbulence. In addition, it may relate to 
shear flow driven instabilities as observed in cloud street formation (Kelly 1994), Taylor-Couette 
flow (Andereck et al. 1986), Poisseuille-Benard convection (Kelly 1977; Yu et al. 1997; Miiller 
et al. 1992), or turbulent bursting in Couette flow (Bottin et al. 1998). 

3. Experiment 

We conducted our experiments in a thin layer of compressed CO2 within a rectangular cell 
of height d = (388 ± 2)/zm and dimensions (100 x 203)e2. The C0 2 gas was 99.99% pure and 
the pressure was (56.5 ± 0.01) bar regulated to ±55 mbar. The mean temperature of the con- 
vection cell was held constant at (28 ± 0.05)°C regulated to ±0.3 mK. For these conditions, the 
Prandtl number was P = 1.1 ±0.04 as determined from a materials property program (de Bruyn 
et al. 1996). The characteristic vertical diffusion time was t v = d 2 / k = (1.532 ± 0.015) sec. 
The experiments were conducted at a fixed inclination angle of 7 = (30.00 ± 0.02)°. We ob- 
tained images of the convection pattern using a digital CCD camera, via the usual shadowgraph 
technique (de Bruyn et al. 1996; Trainoff and Cannell 2002). 

We collected data at 17 equally spaced values of e between 0.04 and 0.22, each reached by 
quasistatic temperature increases from below. At each e, we recorded at least 400 series of 100 
images (at 3 frames per second), with each series separated from others by at least lOOrt, to guar- 
antee statistical independence. These runs are later referred to as short runs. For the lowest value 
of e — 0.04, we recorded up to 600 series to cope with the reduced number of defects present 
in that regime. In addition, we collected data while decreasing the temperature quasistatically 
from e = 0.12 to e = 0.06 to check for hysteresis, which was not observed. Particularly long 
time series were recorded at e = 0.08 (one run) and 0.17 (two runs). These consisted of six-hour 
(1.4 x lO 4 ^) continuous runs of data at 3 frames per second. These runs are later referred to as 
long runs. 

Figure 4a shows a sample Fourier-filtered shadowgraph image of the complete experimental 
cell at e = 0.08. While on first sight the pattern appears to be uniformly composed of UC, 
detailed investigation demonstrates that this is not the case. In Fig. 4b we have superimposed 400 
statistically independent images at the same e (one from each short run). Clearly visible are gray 
regions, such as within the white-framed subregion, where the average looks almost structureless. 
There exist also long-term ordered regions such as the strong black/white roll patches near the 
boundaries. Other spatiotemporally chaotic systems (Ning et al. 1993; Gluckman et al. 1993) 
also exhibit similar behavior due to the presence of boundaries. 

Another method of characterizing the spatiotemporal homogeneity of a "chaotic" pattern is to 
superimpose defect trajectories obtained for statistically independent frames (see Section 5.1 for 
details on defect-detection). An example is shown in Fig. 4c, where we superimpose defect tra- 
jectories for 50 statistically independent short runs. Defects with positive (negative) topological 
charge are white (black). Again, far away from the boundaries (within the marked subregion), 
persistent and homogeneous creation and annihilation of point defects was observed. In contrast, 
other regions were either rarely traversed by defects or were preferentially traversed by defects 
of like charge. In addition, favored spots along the boundaries preferentially seeded defects of 
the same charge, which then followed similar trajectories. Thus, to minimize the impact these 
effects and utilize as ideal a cell as possible, the data presented in this paper were taken only 
from the subregion of dimension 51d x 63d marked by the white box. 
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FIGURE 4. (a) Example convection cell, (b) average of 400 statistically independent images, and (c) defect 
trajectories for fifty statistically independent short runs. White (black) are defects with positive (negative) 
topological charge. The nearly homogeneous subregion (5 Id x 63d) used in all analyses is marked by a 
white box. Uphill is at the left, e = 0.08, P = 1.1 ± 0.04 and 7 = 30°. 

4. Theoretical description 

We base our theoretical description on the standard Oberbeck-Boussinesq equations (OBE), 
slightly generalized to cover the inclination of the fluid layer. As in the experiments, we con- 
centrate on systems with large aspect ratio, i.e., with the lateral dimensions of the cell in the 
x, y plane » its height d in the z-direction. Thus, we adopt the usual idealization of periodic 
horizontal boundary conditions, which are incorporated numerically by switching from position 
space to 2-dimensional Fourier space. 

First, the onset of ILC is obtained by a standard linear stability analysis of the base state. Then, 
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in the nonlinear regime, we investigate the stability of periodic solutions using both weakly non- 
linear and multimode Galerkin approaches. Most important for the interpretation of the (weakly) 
disordered states in the experiments have been fully three-dimensional solutions of the OBE. 

4. 1 . Basic equations 

As shown in Fig. 2, the convection cell is subjected to a temperature difference AT = (T\ — 
T 2 ) > between the hot (bottom) and cold (top) plate and is inclined by an angle 7. In a coor- 
dinate system aligned with the cell, the gravity vector is given by g = — g (zcos7 + y sin 7), 
where g is the acceleration due to gravity. As usual we introduce the Rayleigh number R as a 
nondimensionless measure for AT: 

R ^ C ° S 7 ATd3 (4.1) 

VK 

with a the thermal expansion coefficient, v the kinematic viscosity, k the thermal diffusivity. 
Please note that R depends on 7. 

The OBE become dimensionless if lengths are measured in units of d, time in units of the 
thermal diffusion time t v = d 2 / 'k, temperatures in units of AT/R and velocities in units of 
d/r v . Thus, we arrive at the following nondimensional form of the OBE: 



V 2 u + 9 z + {-Rz + 9) y tan7 - Vtt 



-u+(u.V)u 



(4.2a) 



V 2 6» + i?u-z = — + (uV)0. (4.2b) 

(AT/R)9 denotes the deviation from the linear temperature profile \{T\ + T2) — ATz 
(—1/2 < z < 1/2) of the basic non-convecting state. Terms which can be expressed as gra- 
dients are included in the pressure term \7n. We assume incompressibility of the velocity field, 
u, i.e. V • u = 0. Furthermore 9 and u are required to vanish at z = ±1/2 (rigid boundary 
conditions). 

Even below the onset of convection, the OBE provide a base flow with a cubic mean-flow 
profile TJo(z) (in dimensionless units): 

TT . . i?tan7 / , z\ „ 

U (z) - — g-^ (z 3 - -J y, (4.3) 

as illustrated in Fig. 2. Above the onset of convection, it is useful to split the solenoidal velocity 
field u into three parts using 

u = U (z)+v(a;,y,z;t) + U(z,t). (4.4) 

The term v(x, y, z; t) is derived from a poloidal-toroidal decomposition of u in the form v = 
v po ; + v tor = V x V x (zx) lVx( with the poloidal (and toroidal) velocity potentials \ (and 
Q. This decomposition automatically fulfills V • v = 0. The field v po ;(x, y, z; t) describes the 
periodic convection rolls with their (nearly circular) streamlines and vanishing vertical vorticity 
O z = (V x u) z , while w tor is associated with finite Q z . 

Convection also leads to a modification U(z, t) (see Eq.4.4) of the basic shear flow profile 
Uo(z). The governing equations for U(z, t) = (U x (z, t), U y (z, t), 0), derive from a horizontal 
average of Eq. 4.2a and read as follows: 

(8 2 zz - P^dt) U x + P-^.M = 0, (4.5a) 

(d 2 zz - P^dt) U y + p- 1 d z \^ z ]+9sin 1 - 0, (4.5ft) 

Here, the overline denotes an in-plane average (the q = component in Fourier-space) of the 
corresponding terms. U(z, t) shares the odd z-symmetry with U (z), but with the flow direction 
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reversed. Thus, the baseflow \]q{z), which is potentially unstable due to the inflection point at 
z = 0, is reduced with increasing R. Note that Eqs. 4.2 become 7-independent for LR which 
vary only in the x— direction. Thus the onset of LR and their structure (but not their stability) 
can be deduced immediately from the standard RBC problem without inclination (Clever 1973). 

While w tor vanishes in periodic roll patterns, it is excited by any roll curvature or defects in the 
pattern. In contrast to v po i, the Hagen-Poisseuille-like z— profile of \r tor is even in z, and is often 
referred to as mean-flow. Other fields can thus be efficiently advected by v tor , typically leading 
to further enhancement of some initial perturbation. For Prandtl number of order unity and below, 
the resulting positive feedback is the source of spatiotemporal complexity as observed in defect 
turbulence (Daniels and Bodenschatz 2002a) and spiral defect chaos (Morris et al. 1993). 

4.2. Periodic solutions 

It is convenient to substitute Eq. 4.4 into Eqs. 4.2 and to use a condensed notation for the resulting 
OBE equations: 



The fields x, (, 9, u have been collected into a symbolical vector V. The letters C and C represent 
the resulting linear operators and J\f the quadratic nonlinearities in Eqs. 4.2 and 4.5. 

The first step in any theoretical investigation is the linear stability analysis of the base state 
V = 0. Neglecting TV in Eq. 4.6, we arrive at a linear eigenvalue problem which diagonalizes 
in Fourier space via the ansatz V(x, z, t) = e At e lq x V;i„(q, z) with x = (x, y) and q = (q,p). 
The eigenvalue X(q, R) = a + iuj with the largest real part provides the growth rate a(q, R). 
The condition a(q, R) = describes the neutral stability curve Ro(q) which has a minimum 
at R c = Ro(q c ) for the critical wave vector q c . In the present case, we find uo = 0, so the 
bifurcation is stationary. 

In contrast to isotropic systems, where cr(q) and i?o(q) depend only on the modulus of q, 
the anisotropy causes them to additionally depend on the angle between q and y. Inspection of 
Eq. 4.2 shows that the control parameter R appears not only explicitly in the buoyancy term of 
Eq. 4.2b, but also in the coefficient P _1 i? tan 7 arising from the U contributions (see 4.3) to 
the advection terms (see Eq. 4.4). 

Where buoyancy prevails (7 less than the critical angle 7 C 2(-P)), we find LR at the onset of 
convection, with q c = (q c ,0). The threshold R c = 1708 (defined according to Eq. 4.1) and 
q c = 3.117 are 7-independent. In contrast, for 7 > 7 c2 the instability is shear-flow driven and 
TR with q c = (0,p c ) are selected. This critical behavior at threshold was discussed in Hart 
(1971a,b); Clever (1973); Ruth et al. (1980); Clever and Busse (1977). 

We use Galerkin methods to study ideal periodic convection patterns in the nonlinear regime. 
All fields are expanded in two-dimensional Fourier series with respect to the (x, y)-coordinate 
plane and in suitable test functions in z fulfilling the rigid boundary conditions at z = ±1/2. For 
instance, the component of V reads 



with S n (z) = sin(n7r(,2 + 1/2)). To ensure real-valued 8, we require c„ fe _ m;n = c* k m . n . The 
fields x, C are analogously represented except that we use the Chandrasekhar functions C n (z) 
(Chandrasekhar 1961) instead of the S n (z) for \. 

For LR with q || x all expansion coefficients vanish for m / and conversely for k = 
for TR with q || y. By inserting Eq. 4.7 and the analogous expansions for \, u into Eq. 4.6 
and truncating the series one arrives at a nonlinear algebraic system for the their expansion co- 
efficients. This system of equations is solved using a Newton-Raphson iteration scheme and 



C—V(x,z,t) = CV(x,z,t) +JV(V,V). 



(4.6) 




(4.7) 
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subsequently tested for stability with respect to linear perturbations. Using the standard Floquet 
ansatz, the 9 perturbation reads for instance 



<W(x, z; t) = e a " l e lsx 



cxp [i(kqx + mpy)] ^ Set 



(4.8) 



with s = (s x , s y ). Positive nonlinear growth rates 5?e[cr„(q, s, R, P)] signal secondary instabil- 
ities of the periodic convection patterns. Numerical results agree with all those published previ- 
ously in the literature (e.g. Clever and Busse (1977)). 

This paper investigates the secondary, modulational, wavy (zig-zag, undulation) instability of 
LR at R = R u (-f,P). This instability is initially characterized by transverse undulations along 
the roll axis at large wavelength (s x — and s y — > 0) as first described in (Ruth et al. 1980; 
Clever and Busse 1977). It develops for R> R u slightly above R c , i.e. at a small reduced control 
parameter e u (P, 7) = (R u — R c )/R c . For example, e u — 0.016 for 7 = 30° on the dashed line 
in Fig. 3. 

Weakly nonlinear analysis provides some insight into the underlying physical mechanisms 
of the instability, particularly the final-amplitude state. We utilize a set of coupled, complex 
amplitude equations, which can be expected to work reliably for small e. 

We use the following ansatz for the solution of Eq. 4.6: 

V(x,z,t) = A(t)V ;m (q 1 ,z)e^ x +B(t)V ;m (q 2 ,z) e ^ x -C(t)V ;m (q 3 ,z)e^ x +c.c. + h.o.t. 

(4.9) 

with qi = (q, 0), q 2 = (q,p) and q 3 = (q, —p). The first term describes LR of wavenumber 
q and complex amplitude A(t), while the additional terms provide transverse modulations with 
wavenumbers ±p and complex amplitudes B(t) and C(t). Note that the wavevectors q^ span a 
resonant tetrad obeying 

qi=q2+q3-qi- (4.10) 
We insert Eq. 4.9 into Eq. 4.6 and retain terms up to cubic order in the amplitudes A, B, C. 
After projecting these equations onto the linear eigenvectors Vu n (qj , z), we obtain the following 
coupled amplitude equations, valid for R> R u . 



(4.11a) 
(4.11/7) 
(4.11c) 

The linear growth rates are defined as: <ta(R) = c(qi, R), <Jb(R) — cc(-R) = (7 (q2, R) < 
a A- The nonlinear coupling coefficients bij > and the cross-coupling coefficients pi > which 
characterize a resonant tetrad coupling were calculated using the usual projection techniques. It 
is convenient to characterize the complex amplitudes by their moduli and phases: 

A=\A\exp(i4> A ), B= |S|cx P (z0 B ), C = \C\ exp {icj> c ). (4.12) 

Two phases can be arbitrarily chosen due to translational invariance in the x, y directions and 
the choice <j> A — 0, <f>s — leading to real A, B is convenient. The undulation growth rate 
is obtained by a linear stability analysis with respect to B and C. Solutions with \A\ > and 
\B\ = \C\ > exist only if the condition a u (R) = <T B (R) + (crA(R)+cos((f)c)p2+b2i)/bii > 
holds. <r u (R) is minimal for (f>c = since p 2 > and crosses zero at a value R = R u , which 
agrees with the full Galerkin result. 



dt 


(<TA 




- b\2 


B\ 2 


~b 12 \C\ 2 )A + Pl A*BC, 


dt 




-b 2l \A? 


- b 22 


B\ 2 


- b 23 \C\ 2 ) B + P2 AAC* , 


dt 




-fo 3 i|A| 2 


- b 32 


Bf 


- b 33 \Cf) C + P3 AAB* . 
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FIGURE 5. Stability islands for ordered undulations at P — 1.14 and 7 = 30° in the q, p plane at 

increasing values of e. 



In more physical terms, the undulation solution creates horizontal variations of the mid-plane 
temperature field, which correspond to the experimentally-obtained shadowgraph pictures. It is 
convenient to normalize the temperature component of V;i„(q, z) to unity at the midplane (z = 
0). Thus, we arrive from Eq. 4.9 at the the following representation for the mid-plane temperature 
field of an ideal undulated pattern: 

^ R (x, y) = K [e i<lx (A(q,p) + B(q,p)e^ - C(q,p)e- ipy ) + h.o.t.)] . (4.13) 

with real A, B and C = B and higher order terms neglected. It is illuminating that, slightly 
above the instability, ipR as presented in Eq. 4.13 can be considered as a phase modulation of 
the LR pattern. In fact, a modulated LR pattern 2 A cos(qx + <j>{y)) with <p(y) = ^ sin(py) and 
j<1 corresponds to Eq. 4. 13 for the resonant case B = C in a leading order expansion with 
respect to B/A. 

The amplitudes A and B = C obtained from Eqs. 4.1 1 and the resulting approximate solution 
V according to Eq. 4.9 are convenient for generating the full nonlinear Newton iteration scheme. 
In this way, we construct fully three-dimensional ordered undulation (OU) solutions of the OBE 
of the type described in Eq. 4.7 for R> R u . The subsequent linear stability analysis (see Eq. 4.8) 
yields the stability regions of the OU. 

Figure 5 shows the representative results for P = 1.14 and 9 = 30°. The OU at a given e are 
stable within certain regions (islands) in the (q, p)-space. These islands shrink with decreasing e, 
becoming arbitrarily small as e approaches e u . In this limit, the island is located at q w q c and 
p w and its area approaches zero. The islands reach their maximal size for 0.1 < e < 0.2, 
before they shrink again as the transition to crawling rolls is approached. 
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FIGURE 6. Simulation at 7 = 30 and e = 0.025 (R = 1750). (a) Stationary OU after 10 5 iterations 
(2000^), starting from random initial conditions, (b) Fourier spectrum of (a): pattern locked perfectly into 
the resonant tetrad described by the amplitude equations. 

The limits of the islands in Fig. 5 are given by amplitude instabilities (s = in Eq. 4.8) at the 
low-p boundary of the islands and long-wavelength modulational instabilities (\s\ <C q c ) at the 
high-p boundary. In the former case, the amplitude A is virtually unaffected, whereas the other 
amplitudes in Eq. 4.13 are perturbed as follows: B ^ B + 5, C ^ C — S, with \6\ <C B. Due to 
the asymmetry B^C, a non-vanishing component U x develops, i.e. a mean-flow which likely 
drives the defect-turbulent state. By analogy to the aligning effect of the basic profile Uo on the 
LR, U x will tend to turn the roll axis as well. 

In contrast to the amplitude instabilities, the modulational instabilities (which involve finite 
Floquet vectors s) will change the basic periodicity of the undulation solutions. The main con- 
sequence of such a perturbation involves a mean-flow (along x) contribution to v tor , which 
presumably induces defect chaos as well. 

4.3. Simulations 

Since the Galerkin method is restricted to the analysis of stationary periodic states and their 
stability, direct numerical simulations of the OBE is needed to investigate the time evolution 
of complex patterns. For that purpose, a previously-developed code (Decker et al. 1994; Pesch 
1996) was generalized to cover inclination (Brausch 2001). This code represents all fields by an 
appropriate Galerkin-ansatz like Eq. 4.7 and treats the (x, y) -dependence in Fourier space by a 
pseudo-spectral technique. By using a small (< 10) number of modes in the z direction, reliable 
simulations of large aspect ratio systems can be conducted on common computer clusters in a 
reasonable time. In the (x, y)-plane, our simulations covered areas up to 40c? x 75d, where we 
used up to 384 Fourier modes in either direction. Periodic boundary conditions have been used 
except in a few cases, where we introduced subcritical e ramps to suppress convection at the 
sidewalls to approximate experimental conditions. We typically found the same bulk behavior in 
both cases for simulation domains of size 40c? x 40c? or greater. 

The time integration in the simulations was started either from random initial conditions (to 
scan the manifold of nonlinearly selected solutions) or from previously calculated Galerkin solu- 
tions with superimposed noise to check for stability. A representative example of stationary OU 
is shown in Fig. 6; the system chose wavenumbers q = 1.05q c = 3.26 and p — 0.16g c = 0.5, 
which is inside the corresponding stability island in Fig. 5. The resulting ratio B/A = 0.59 
according to Eq. 4.13 agrees with Galerkin results. In fact, OU were always selected for small 
e < 0.03 when starting from random initial conditions. At larger e OU could still be obtained 
by starting from a stable periodic Galerkin OU with a small amount of superimposed noise with 
a signal-to-noise ratio of 10~ 4 . This is demonstrated in Fig. 7a. Here we started at e = 0.1 and 
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FIGURE 7. Convection patterns at e = 0.10 and 7 = 30°. (a) Numerical simulation started from Galerkin 
solution with superimposed Gaussian noise, showing perfect OU after reaching steady state, (b) Numerical 
simulation started from Gaussian noise, with subcritical e-ramp simulating lateral boundaries. OU are vis- 
ible in the upper left and UC in the lower right, (c) Fourier-filtered shadowgraph image of homogeneous 
subregion of experimental cell. 




FIGURE 8. Five views of a demodulated experimental convection pattern showing undulations and defects, 
(a) Real part i\>b. of a Fourier-filtered shadowgraph image at e = 0.10 with negative (left) and positive 
(right) defects. Uphill direction is at the left side of the page, (b) Phase field <j>: black is <j> = and white 
is <t> — 2ty. The arrow indicates an integration loop around a defect as given in Eq. 5.1. (c) Zero crossings 
of ipn and ipi\ white is tpa = and black is ipi — 0. (d) Modulus of image, \ip\. (e) Power spectrum of 
ipR. (q,p) are wavenumbers associated with (x, y) directions. Shaded region of Fourier space (q < 0) is 
discarded during reconstruction of the complex field. 

q = 1.05q c = 3.28, p = 0.31q c = 1.02 inside the corresponding stability region in Fig. 5 and 
recovered the ideal Galerkin solution with B/A = 0.94. However, the basin of attraction of the 
OU shrinks with increasing e and the simulations starting from random initial conditions typi- 
cally settle into UC, as for example shown in Fig. 7b. The ensuing bistability between stationary 
perfectly ordered solutions and dynamic weakly defect turbulent at larger e is the main theme of 
this paper. 



5. Analysis of complex undulation patterns 

In this section we discuss several approaches to characterizing STC in general and the onset of 
chaotic undulations in particular. Section 5.1 explains how we identify all topological defects in 
a shadowgraph image, and generate a time series of the total number of defects, N(t), from the 
consecutive images of an experimental run. Section 5.2 deals with the pattern entropy S(t) and 
the transverse correlation length Finally, in Sec. 5.3 we discuss the statistical properties of 
the amplitudes A, B, C, defined in Eq. 4.9. 
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5.1. Identification of defects 

A shadowgraph (or simulation) image can be understood as the real part ipR of a complex field 
tjj(x, y) — \ip(x, y)\ exp[i(j)(x, y)] with modulus \ip\ and phase cf). To construct ip(x, y), we de- 
modulate the image as illustrated in Fig. 8: after a 2d Fourier transformation, half of the modes 
in the Fourier plane are set to zero (see Fig. 8e). A subsequent inverse Fourier transformation 
recovers tp(x, y). 

As already discussed above, the characteristic feature of a UC pattern is undulating stripes 
containing topological point defects. Defects locally change the spacing of rolls and their orien- 
tation (Bodenschatz et al. 1988). Complex analysis requires that defects be located at a zero of 
ip(x, y): i.e. at the points where the lines of zero real part, iPr(x, y) = 0, and zero imaginary part, 
ipi(x, y) = 0, cross (Rasenat et al. 1990). (See Fig. 8bc for an illustration.) As a second criterion, 
the topological charge (or "winding number"), n, defined by the circulation of the phase gradient 
about a defect 

<j)V<t>-ds = ±n2ir (5.1) 

becomes nonzero; only point defects with \n\ = 1 are found in UC. 

To locate the zero crossings efficiently (Egolf 1999), we formed a set of four binary fields for 
4>R{i,j) at each pixel ((i,j) <G Z) of the plane. 

fa^\i,j)={ 1 ' ^ + A ^ + A -?)>° (5.2) 
0, otherwise 



where (Ai, Aj) <G {(0, 0), (0, 1), (1, 0), (1, 1)} denote the neighboring pixels of (i, j). The four 
binary fields ip R are then logically combined at each point to produce a new binary field 

^r = (fa © fa 1 ) v (fa e fa ) v (fa © fa 1 ) (5.3) 

where © is the XOR (exclusive OR) logical operator and V is the OR logical operator. Anal- 
ogously, we define the new binary field ipj{i,j) as well. All points (i, j) where ip R (i,j) = 
?Pi(i,j) = 1 serve as candidates for a defect. To ensure that such a point is indeed a defect, 
the contour integral in Eq. 5.1 is evaluated by traversing its eight nearest neighbors. A defect 
candidate is accepted if the contour integral yields ±2ir, assigning the appropriate sign to the de- 
fect based on the sign of the integral. After locating all defects, we eliminated double-detections. 
Counting these zeros of ij)(x, y) for a time series of frames yields N(t). 

5.2. Global stochastic properties of undulation chaos 

The following section continues to analyze the representative case of inclination angle 7 = 30°. 
A prominent feature of the experiments is that at higher e > 0.10 the competing attractors - OU 
with few defects and UC with many as shown in Fig. 1 - are alternately visited in time, while the 
system remains persistently chaotic at lower e. A similar bistability was also found in the case 
of (isotropic) RBC: Cakmur et al. (1997) observed competition between ideal rolls (IR) and the 
spatiotemporally chaotic state of spiral defect chaos (SDC). In that case, however, a transition 
from SDC to IR was observed by decreasing e. There, the analysis of the spectral entropy S(t) 
(Neufeld and Friedrich 1994) was illuminating, and the concept is applied here as well. S(t) is 
defined as 

S(t) = -(P(q,p,t)lnP(q,p,t)) (5.4) 

where P(q,p, t) is the normalized spectral distribution function that describes the power in the 
mode with wavevector (q,p) at time t and the average (•) is performed over the (q,p) plane. 

The function S(t) provides a measure for order in a pattern: we have S — In 6 for OU de- 
scribed by the six modes of Eq. 4.13 while a disordered pattern has more modes excited and 
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FIGURE 9. Spatial autocorrelations (M(r)M(r — ro)} of the modulus M — \ip\ of the complex demod- 
ulation functions xf> associated with the images shown in Fig. 1. (a) UC and (b) OU. Images are rotated by 
6° clockwise with respect to the coordinate system shown in Fig. 1, due to a slight misalignment of the 
experimental cell. 




6000 7000 




FIGURE 10. Time trace of (a,d) spectral entropy S, (b,e) number of defects N, and (c,f) transverse cor- 
relation length £ at e = 0.08 (a,b,c) and e = 0.17 (d,e,f), extracted from shadowgraph images. Data was 
sampled every 10 seconds (6.5r„). 



S > In 6. Conceptionally for anisotropic patterns, like undulations, the pattern entropy is even 
better suited to characterize disordered pattern than in isotropic systems. For instance a target 
pattern is represented by a ring Fourier space and therefore has a deceptive high value S, though 
it is a well ordered state. 

Further characterization of the relative order of the pattern is possible using the correlation 
length £ in the transverse (y) direction obtained from the half-width at half-maximum of the spa- 
tial autocorrelation functions, as demonstrated in Fig. 9. In image (a), the pattern is disordered, 
the spatial autocorrelation function is sharply peaked at the origin, and £ only a few d. In image 
(b), the pattern is highly ordered and shows correlations that fall off with a correlation length 
approaching the length of the cell. 

Fig. 10 shows example time traces of N(t), S(t) and for e = 0.08 in the UC regime 
and also for e = 0.17, well above the transition to competition between UC and OU. The dashed 
vertical lines for the case e = 0.17 mark the times of the images shown in Fig. 1 , which exemplify 
the cases of high S (disordered UC) and low S (ordered OU), respectively. For each time series of 
the quantities X = S(t), N(t),£(t) we obtain the time average (X) and the standard deviations 

e x = v^ 2 ) - {xy. 

For e = 0.08 the quantities S(t), N(t), and £(t) show smaller fluctuations than at e = 0.17. 
In fact, at higher e the system visits alternately OU and UC states, and we observed a num- 
ber of instances where the system remained in the OU state for more than 500 r„. Because 
ordered patterns have few defects and chaotic patterns have many, one expects these measures 
to be correlated with each other. Indeed, the number of defects N exhibits behavior mirroring 
that of S(t) while £(t) is negatively correlated. Figure 11 shows temporal cross-correlations 



14 



K. E. Daniels, O. Brausch, W. Pesch, and E. Bodenschatz 





-1000 -500 500 1000 

At W 



-1000 -500 500 1000 




-1000 -500 500 1000 



FIGURE 11. Cross-correlations C(S, N), C(7V,£), C(£, S) of spectral entropy S(t), number of defects 
N(t), and correlation length from left to right. 
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FIGURE 12. (a) Mean correlation length (£) and (b) its standard deviation as a function of e. (c) Mean 
spectral entropy (S), and (d) its standard deviation Es as a function of e. Solid lines oc \ft — e u are plotted 
as guides to the eye. Triangles are for short runs in which temperature steps between them were either 
increasing (A) or decreasing (T). 



C x ,y(At) = ((X(t + At) - (X))(Y(t) - for X, Y e S(t),N(t), and In 

all cases, the data at e = 0.17 is more strongly correlated. 

The transition to increasing order, i.e. the appearance of long-living OU states for e w 0.1 
is a peculiar one, in that the system does not remain in the new, ordered state but intermittently 
returns to the chaotic state. To characterize this transition we measured the quantities (S), (£), Eg 
and at each e, shown in Fig. 12. All four plots exhibit sudden changes near a transition point 
e t s=s 0.1. For instance (£)(e) increases sharply there, indicating that the size of ordered regions 
is strongly growing for e > e t . Both £s(e) and S^(e) show that the fluctuations increase in 
magnitude as well when e > e t . No hysteresis was observed across this transition. Note that this 
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transition also corresponds to qualitative changes in the defect dynamics as previously described 
in Daniels and Bodenschatz (2002a). 



The analyses presented above are not specific to undulation chaos. Thus, we have used an ad- 
ditional method, to characterize slightly disordered undulation patterns. The idea is to use the 
ansatz in Eq. 4.13 locally, i.e. with real space- and time-dependent quantities (q,p, A, B, C). 
The first step is to extract from our experimental images the complex function t(j(x,y) as ex- 
plained at the beginning of the present section (see Fig.. 8d). Then q(x,y) is obtained by the 
local wavenumber method described in Egolf et al. (1998), by calculating 



Near zeros of ip(x, y) , for instance at defects, the wavevector field q(x, y) field is ill-defined 
and yields unphysically large values, which can safely be identified and neglected. To determine 
the local undulation wavenumber p(x, y) we use the squared modulus of ip, which can be locally 
well-parametrized by the following analytical expression (see Eq. 4.13): 



Consequently the local wavenumber analysis now applied to \ip\ 2 yields p(x,y). For weakly 
disordered OU the periodicity of \t/j\ 2 is in fact governed by 2p. For B s=a C, Eq. 5.6 reduces to 



Fig. 13 shows sample wavenumber distributions of (q,p) for the experimental convection pat- 
terns cells together with the stability islands from Fig. 5. The mean values of q,p (the small 
circle) are in the stable regime. We previously observed the mean value of p to be proportional 
to ^Je — e u , where e u agrees with the Galerkin prediction for the onset of OU (Daniels and Bo- 
denschatz 2002a). The mean roll wavenumber q does not vary significantly with e. Inspection of 
Fig. 13a shows that the distribution of the q, p values overlaps with the stability regions in a wide 
range, and the most probable value (circle) is within the stable regime. Near defects, the undu- 
lation patterns are strongly deformed and local wavenumbers can be pushed outside the stability 
regions. 

Further insight arises from distributions of the amplitudes A, B, and C. For an ideal undulation 
pattern characterized by the wavenumbers q, p we also know the amplitudes A and B = C ac- 
cording Eq. 4.1 1. As demonstrated in Fig. 13b the ratio C/A varies significantly with wavenum- 
ber. Thus, it is not surprising that a disordered pattern with a fairly wide range of q, p values (as 
shown in Fig. 13) is characterized by a broad distribution of the amplitudes A, B ^ C as well. 

We determined A, B, and C from the local maxima and minima of the experimental \ip\ 2 . Due 
to the presence of defects it is not feasible to keep track of the phases of the cosine functions in 
Eq. 5.6. Since y — > y = y + n/p results in B <-> C we have chosen the convention B < C in 
all cases. Our analysis is independent of q(x, y) and p(x, y) and is not sensitive to the tilt in the 
undulations. As shown in Fig. 14, the median value of | ^ | is observed to be low (0.2) near onset 
and rise with e. For small e the peak in 1 1| | is near the ideal OU value B/C = 1. With increasing 
e the intermittently-visited UC states are responsible for an increasing asymmetry B ^ C, which 
is then reflected in the average ratio in Fig. 14 as well. 

The averages of the ratios C/A and B /C shown in Fig. 14 are obtained from a sequence of the 
corresponding distribution functions examples of which are shown in Fig. 15. We observe that 
the distributions become slightly more narrow at the higher e, which might indicate the increased 
prevalence of OU. 



5.3. Local amplitude analysis 




(5.5) 



|Vf = A 2 + B 2 + C 2 — 2BC cos (2py) + 2A(B - C) cos (py). 



(5.6) 



\%P\ 2 = A 2 + 2C 2 - 2C 2 cos (2py). 



(5.7) 
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FIGURE 13. (a) Stability islands at e = 0.1 from Fig. 5, showing sample wavenumber distributions for each 
of the images in Fig. 7. Peaks of the distributions are given by symbols; lines are the half-maximum of the 
distributions showing UC. (b) Corresponding theoretical amplitude ratios |§| for q = 3.1, 3.3, 3.5 within 
the stable region. 




FIGURE 14. Ratios | ^ | and 1 1| | as a function of e, measured at the peak value of the histogram. Symbols 
represent averages from 500 experimental images, for which quasistatic temperature steps are increasing 
for A and decreasing for T. 




FIGURE 1 5 . Histograms of j ^ j and j || | at two values of e. Averaged experimental points (solid) are taken 
over 500 images, with comparison to single images (dashed, dotted) as specified in legend for e = 0.17. 
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One feature not relevant for the present analysis is a weak drift of the experimental patterns 
down the cell (Daniels and Bodenschatz 2002a). This is also the reason that undulations are 
mostly invisible in the time-averaged image shown in Fig. 4b. We associate the drifting behav- 
ior with non-Boussinesq effects, which combined with the effect of the basic shear render the 
growthrates of (Jb,&c of the oblique modes in Eqs .4. life and 4. 11c complex. Between e = 0.08 
and e = 0.17 a drift speed of about O.Q5d/r v has been measured, which is is substantially smaller 
than the value d/r v w 1 observed above -f c2 . The difference may be partially accounted for by 
the strength of the shear flow, which is greater for steeper inclinations. 

6. Discussion 

The central topic in this paper is the competition between ordered (OU) undulations and un- 
dulation chaos (UC). The Galerkin analysis of the OBE demonstrates the existence of linearly 
stable undulation patterns above the wavy instability up to fairly large e. As shown in Fig. 5 
the wavenumbers q,p, which characterize the three dominant modes A, B, C (Eq. 4.13) cover 
a finite region in q-space, which is maximal near e = 0.1. The linear stability of OU is con- 
firmed by direct simulations of the OBE when starting from a Galerkin solution with noise (see 
Fig. 7a). However, the basin of attraction of the OU shrinks with increasing e. At small e starting 
from random initial conditions the system selects OU (see Fig. 7a), albeit after a long transient. 
In contrast, at larger e the system seems to remain in a UC state that competes with regions of 
OU (see Fig. 7b). Within the spatial and temporal limitations of our simulation we are unable to 
determine whether we observe UC as a long-lived transient or as a truly asymptotically selected 
state. 

In the experiment we do not observe ideal OU in any regime. At small e, for instance at 
e = 0.08, the number of defects N in our sample area (the box in Fig. 4) varies between 3 and 
18; the average is about N — 10 with fluctuations between these extrema (see Fig. 10b). For 
larger e the situation is different (see Fig. lOd). We find an intermittent switching between states 
with very few defects (~ 4 in Fig. lb) and others with many (up to N > 25) as in in Fig. la. 
On the other hand, the system can remain for long times in either of these states, switching 
intermittently between them. 

We interpret this scenario as reflecting the competition between OU and UC states. Defects 
distort the patterns and can easily lead to the regions of the cells where part of the the wavenum- 
ber distribution lies outside the stability regions. The creation of defects and their dynamics is a 
well-known mechanism allowing the system to return to the stable regime. In the optimal case, 
this is achieved when oppositely charged defects annihilate each other or move to the boundaries. 
This process takes time and is susceptible to imperfections in the cell (particularly the bound- 
aries), but there are also cases (such as SDC, (Bodenschatz et al. 2000)) where the system never 
settles down to the stable attractor but switches intermittently back and forth between almost 
perfect ordered and strongly disordered states as in our case. 

7. Conclusion 

In this paper we have studied an anisotropic pattern forming system which shows spatiotem- 
poral complexity. In contrast to the similar and previously studied scenario of the electrohy- 
drodynamic instability in nematic liquid crystals (Rasenat et al. 1990), the basic equations are 
simple enough to be amenable to a precise theoretical analysis. We demonstrate the crucial im- 
portance of topological defects in the competition between chaotic and well-ordered patterns. 
While numerical simulations demonstrate stable ordered undulations at all e above the secondary 
instability, such a state is only intermittently accessible in experiments. Interestingly, this ordered 
state occurs only at higher driving, where the stability region for the undulations is largest. 
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